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Abstract 

Wc analyze the (1+1) dimensional QCD (QCD2) at finite density to consider a 
number of qualitative issues: confinement in dense quark matter, the chiral symme- 
try breaking near the Fermi surface, the relation between chiral spirals and quark 
number density, and a possibility of the spontaneous flavor symmetry breaking. We 
argue that while the free energy is dominated by perturbative quarks, confined ex- 
citations at zero density can persist up to high density. So quark matter in QCD2 
is an example of Quarkyonic matter. The non-Abelian bosonization and associated 
charge- flavor-color separation are mainly used in order to clarify basic structures of 
QCDo at finite donsit>^ 



1 Introduction 

The confinement is a distinct phenomenon in Quantum Chromodynamics 
(QCD). The confinement reveals its property in the physical spectra saturated 
by the color singlet states. Many features of spectra are qualitatively explained 
by the picture of the squeezed colored flux connecting colored objects. 

Recently, roles of confinement in cold, dense quark matter have been addressed 
by McLerran and Pisarski [1] . They suggested that excitations can remain con- 
fined even after quarks are released from baryons to form quark matter. Such 
quark matter is named Quarkyonic matter [1,2,3,4], which is distinguished 
from the conventional quark matter with deconfined excitations. The under- 
standing of excitations modes is important to consider the phase structure, 
transport phenomena, etc. 

Before suggestions in Ref.[l], confining effects in quark matter has not been 
taken into account seriously. Presumably one of the reasons would be that 
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quarks seem to excite individually after being released from baryons. In ad- 
dition, asymptotic freedom apparently implies that confining forces become 
irrelevant when typical distance among quarks is very small. 

In this paper, we are going to offer one counter example to such reasonings. 
While we expect that conventional reasonings work for most of regions in the 
quark Fermi sea, they do not capture entire physics near the Fermi surface. 
To illustrate the points, we analyze QCD in (1+1) dimensions (QCD2) as a 
theory with confining forces and asymptotic freedom [5,6,7,8,9,10]. We will 
conclude that cold, dense quark matter in QCD2 is an example of Quarkyonic 
matter. 

Our main statement is that the mechanism to suppress the confining gluons 
at finite density is neither asymptotic freedom nor the percolation of decon- 
fined bags, but the medium induced color screening. Since the color charge of 
the system is zero in average, the strength of the screening is determined by 
phase space allowed for colored quarks to virtually fluctuate near the Fermi 
surface. In spatial dimensions, d, larger than one, such phase space increases as 
~ /i'^~^AQCD, so the strength of the screening enhances compared to vacuum 
case as ~ (yu/AqcD)'^"^ (/^, Aqcd are quark chemical potential and nonpertur- 
bative scale, respectively). So deconfinement should take place when quark 
fluctuations become comparable to those of gluons. 



where we have multiplied color factors for quarks, N^. and for gluons, N^. In 
our world, (i = 3, so the gluon sector is modified at /i ~ A^^y^AqcD, as given in 
Ref. [1]. This is parametrically larger than the scale, /x ~ Aqcd? where nuclear 
matter is generated and shortly later is turned into quark matter [1] . In spatial 
one dimension, the phase space at finite density is the same as the vacuum 
case, so screening effects are not enhanced at finite density. It means that if 
the excitations are confined at zero densitjQ they can survive even at very 
high density. We will see that this is indeed the case in our QCD2 studies. 

The second purpose of this paper is to examine mechanisms of the chiral sym- 
metry breaking near the Fermi surface (modulo strong infrared phase fluctu- 
ations peculiar to (1+1) dimensions [12,13,14]), which recently attract much 
attentions in (3+1) dimensional quark matter in the context of inhomoge- 
neous chiral condensates [3,4,15,16,17,18,19]. We also study the relationship 

^ The analyses in this paper do not answer whether excitations are confined at 
zero density. In the large limit, the theory is confining, at least in the presence 
of arbitary small current quark mass [5]. On the other hand, in the large flavor limit 
such that Nc/N{ — )• 0, it is argued that the color sector takes a similar form as the 
Schwinger model — Higgs phenomenon happens due to the screening [11]. 




d-l 




(1) 



2 



among a scale generated by the colored sector, cliiral symmetry and quark 
density. In (1+1) dimensions, their relation can be seen at the operator level 
[20,21,22,23,24,25]. In higher dimensions, although such operator relations do 
not strictly hold, at high density the curvature of the Fermi surface is so small 
that there are certain circumstances in which low dimensional picture of the 
Fermi surface is usefuQ The detailed understanding of QCD2 would provide 
a useful way of thinking for these situations. 

One of the most enlightening approaches to discuss the above issues in QCD2 
is the non-Abelian bosonization [21]. The use of it is motivated by at least two 
reasons. 

First, the bosonized form allows us direct access to the colorless objects such 
as quark number density or chiral density. A number of conclusions can be de- 
rived without explicitly treating the colored objects and confining interactions 
among them, which are not easy to deal with. This is not the case in usual 
fermionic expressions in which we have to built up color singlet quantities 
from single quark propagator. The trouble is that even if we investigate sin- 
gle quark propagator very precisely, it is not stable against residual confining 
interactions [6] . So such information is not always converted into the physical 
results until we construct colorless objects for which residual interactions are 
under control. 

The second reason comes from utility of the charge-color separation, analogous 
to the spin-charge separation in condensed matter system [23,24]. In chiral 
hmit, the quark number and colored sectors decouple as [8,25] 

'S'fermion ^ [V'] + 'S'color[^]5 (2) 

SO that bosonized fields responsible for U{1) quark number (ip) and color 
densities (h) can be treated separately. It means that the color sector remains 
the same at any quark density, since quark chemical potential couples only to 
quark number density, not to color density. Hence confined excitations at zero 
density persist up to arbitary high density. 

The charge-color separation is also useful to get insights of the chiral symmetry 
near the Fermi surface. In chiral limit, the chiral condensate in the bosonized 
form has the following factorized structure. 



^ At high denisty, transverse dynamics along the Fermi surface provides negligible 
contributions, ~ to the energy spectra. This can be used to factorize integral 

equations such as the Schwinger-Dyson equation near the Fermi surface by integrat- 
ing out transverse dynamics, instead of just ignoring it. The resulting equation is 
dimensionally reduced one. Such treatments are done in [3,4,15,16,17]. 
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{i^Li^R) = {iih){e'^). 



(3) 



Since (tr/i) is unaffected by quark density, density effects appear only through 
the exponent (e**^). The color sector serves a massive scale for the amplitude 
of the chiral condensate, while the quark number sector provides a phase 
rotation. Through the operator relations, we will see single baryon accompa- 
nies single chiral spiral. Then it follows that at nearly uniform quark density, 
largely overlapped baryons, which are merged with other baryons, induce a 
lot of chiral spirals with period ~ 1/2/i. This period reflects that the chiral 
condensates are made of co-moving particles and holes near the Fermi surface, 
Pf ~ yU- Similar results have been found in other approaches [26]. We will ar- 
gue that chiral symmetry breaking effects are relevant near the Fermi surface, 
while not in most regions of the Fermi sea. 

This paper is organized as follows. In Sec. 2, we give general remarks on 
the Fermi sea in Quarkyonic matter. Then QCD2 is briefly discussed in the 
fermionic language. We will also explain why we think of the quark Fermi 
sea, instead of the baryonic Fermi sea. In Sec. 3, the basics of the non-Abelian 
bosonization are quickly reviewed. In Sec. 4, one baryon, two baryons, and finite 
baryon density are pedagogically discussed. The relationship between chiral 
spirals and baryons is shown. In Sec. 5, we argue the extension to two flavor 
case. In Sec. 6, terms which enhance or suppress the formaion of chiral spirals 
are classified. Sec. 7 is devoted to the summary and discussions. Throughout 
this paper, we fix the value of NcqI when we change A'^c from three. 



2 General Remarks on the quark Fermi sea 



Before performing explicit calculations, we argue qualitative aspects of the 
Fermi sea in Quarkyonic matter. To emphasize conceptual points, we first con- 
sider very high quark density where the inter quark distance is much smaller 
than a typical distance scale of confinement. Baryons already overlap one an- 
other so that quarks need not belong to particular baryons. The lower density 
will be discussed later in order to emphasize how nuclear and Quarkyonic 
matter are related. 

Here we should specify our terminology used in this paper. The confinement 
and deconfinement will be classified by excitation modes. We think that this 
will give more precise classifications than those based on releasing of quarks 
from baryons. For example, the former distinguishes whether correlation func- 
tions for glueball opearators contain deconfined mutli-gluon spectra or notp] 



We have to admit that practically it is not easy to make such a distinction in 
a soUd way. For instance, in chiral limit, the p spectral function at zero density 



4 




Fig. 1. The scattering of quarks in the region (I), |p| < ^ — Aqcd- (a) Smah momen- 
tum exchange processes, which are forbidden by Pauh blocking, (b) Hard momentum 
exchange processes, which transfer quarks from the region (I) to its outside. (Dotted 
hues separate the regions (I) and (II).) 

What are relationships between asymptotic freedom and quark density? When 
the inter quark distance is small, typical interactions are supposed to be weak 
because of asymptotic freedom. This logic is frequently used to justify weak 
coupling treatments of high density QCD as well as to guarantee the picture 
of the degenerated quark Fermi sea. 

For more precise statements, we divide regions of the Fermi sea into (I) an 
inner region {\p\ < jj, - Aqcd) and (II) a surface region {\p\ > /i - Aqcd) (see 
Fig. 1(a)). Let us clarify properties of interacting quarks in each region. 

In the region (I), Pauli principle prevents quarks from being scattered by 
small momentum exchanges (Fig.l.(a)). Thus quarks in the region (I) have 
little chance to feel nonperturbative effects. Allowed processes are hard scat- 
terings which transfer quarks to the outside of the Fermi sea (Fig.l.(b)). The 
latter may be treated within weak coupling methods. Therefore quarks in the 
region (I) can be described by quasi-particles affected by perturbative self- 
interactions and many body effects. 

The above argument is incomplete, since there are also contributions from 
zero momentum transfer processes, which are allowed inside of the Fermi sea. 
Such processes should be sensitive to the gluon propagator at zero momentum, 
and would be problematic. If we compute the total free energy as a sum of 
uncorrelated single quark contributions, zero momentum transfer diagrams of 
quark loops would have the IR divergence from p = component of the gluon 



includes the cut of multi-pions starting from zero invariant mass. On the other hand, 
we know that the height of multi-pion spectra of 0{N^) is much lower than that 
expected from qq continuum of 0{Nc). Unfortunately this classification becomes 
again ambiguous in the quark-hadron duality region at high energy [27]. 
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a) ^ b) 





Fig. 2. The scattering of quarks in the region (II), |p| > — Aqcd- (a) Small 
momentum exchange processes as sources of nonperturbative phenomena, (b) Hard 
momentum exchange processes. Both processes are allowed. 



propagate^ 

There is a case study for the confining propagator in (3+1) dimensions. It 
indicates that zero momentum contributions are canceled out after consistent 
treatments of the quark self-energy and confining interactions inside of color 
singlet object^J Intuitively, this cancellation may reflect that the gluon prop- 
agator at p = 0, which corresponds to an infinitely long string configuratioq ^ [ 
can not be relevant inside of color singlet objects. Although we are not aware 
of whether this cancellation generically holds in any confining models, we will 
rely on this qualitative interpretation for the nonperturbative part, skipping 
the problem of zero momentum transfer processes. 

Now let us consider the region (II). There soft momentum exchanges (Fig.2.(a)) 
are allowed in the region (II). In addition to perturbative effects (Fig.2.(b)), 
nonperturbative effects on quarks also operate here, and properties of quarks 
near the Fermi surface may be strongly modified from those of quasi-free 
fermions. 

With these pictures for quarks, let us first examine bulk thermodynamic quan- 
tities to which all quarks in the Fermi sea contribute. A schematic picture is 
given in Fig. 3, taking the pressure as an example. At high density, they are 



^ A similar problem has been discussed in case of an electron gas, and usually one 
handles this by including the Debye screening mass to kill the infrared divergences 
from photon propagators. 

^ The simplest way to show the cancellation is to just use the principal value reg- 
ulator for the force for which the propagator takes zero value at p = [6]. In 
Ref.[28], usual IR cutoff scheme is used to argue the cancellation. While two schemes 
provide different quark propagators, they give the same spectra at the level of the 
Bethe-Salpeter equation. Thus the detail of p = is irrelevant as far as consistent 
calculations are performed. 

6 In usual IR cutoff scheme, - f dp T^x^e^P'^^ ^g-'^fi^ + ar + 

0{iJLRr'^), where hr is infinitesimal. The first term is independent of p, and expresses 
a string with infinite length. 
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■(II):- 



■(I):- 



Number of d.o.f. : Large 
Quark properties: Perturbative 



Quark properties: Non-perturbative 



Number of d.o.f. : Small 



Fig. 3. The pressure as an example of bulk quantities. The dominant contribution 
(~ comes from the region (I) including a large number of (perturbative) quarks. 

well saturated by contributions from the region (I), simply because a num- 
ber of quarks is much larger than those in the region (II). The region (I) 
gives contribution of ~ /i^ including perturbative corrections, while we get 
small (nonperturbative) contribution of ~ /i^Ag^^, from the region (II). In 
this sense, the picture of perturbative quarks should work in computations of 
bulk quantities, and should reproduce perturbative results [29]. 

Again, this does not means that the pressure in Quarkyonic matter is built up 
from a sum of non-interacting single quark contributions in a literally sense. 
Interacting diagrams must be included, to avoid the problem of infinitely long 
strings generated from quarks. Our claim is that once the nonperturbative, 
infrared problem is handled, quarks inside of the Fermi sea should look like 
perturbative quarks. 

While the bulk quantities should be well-approximated by perturbative con- 
tributions, physics near the Fermi surface is sensitive to the excitation modes 
and nonperturbative effects. Examples are the phase structures, transport 
phenomena, etc. The issues of confinement and chiral symmetry breaking are 
nontrivial here. 

To examine these pictures, let us take one flavor masselss QCD2 as an example. 
The action of QCD2 with axial gauge fixing is 



where x = {t,z), = iptA'J^'ip, and is usual color matrix in the founda- 
mental representation. Its normalization is tr[tA^_B] = ^ab/'^- The interactions 
between colored currents are confining, and instantaneous in axial gauge, 
D^^{x — y) = S^^6fj_o6^o\x — y\. The gluon propagator is not dynamical in 
(1-1-1) dimensions, so we have eliminated the gluon fields using the equation 



S 
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of motion. 

In the presence of the Fermi sea, it is convenient to redefine quark fields 
depending on their moving directions -\-z or — 

^±(t,^)-e±^'^>;(t,^) . (5) 

Its shorthand notation is ■?/' = e^^'^^^ip', since eigenvalues of 75 = 7o7z coincide 
with moving directions in (1+1) dimensions, 

tpi^tp = ipU{do + dz)ip_ + ip\.i{do - 8^)1^+. (6) 

So Eq.(5) just means shifts of momenta to measure them from the Fermi 
surface, ±//. With these new fields, we obtain 

V;(x)(i^ + /.7°)V'(^) ^ i^'{x)i0i;'(x) , J^{x) ^ J"^,{x), (7) 
and the Lagrangian becomes that of the zero density. 

Since Lagrangian is the same as that in vacuum, properties of excitations 
such as energies are also unchanged except their shifted momenta. Therefore 
in (1+1) dimensions, if models have confined excitations in vacuum, they also 
do at finite density. This is true no matter how quark density is high and 
inter quark distances are short. The value of A'^c is not essential in the present 
discussion. An overlap of baryons and deconfinement do not have one to one 
correspondence. 

In coordinate space, this situation may be difficult to imagine if we view quarks 
as point-like particles as classical mechanics. Perhaps it is easier to understand 
if we see the distribution of colors by focusing on the wave properties of quarks. 
Closely packed quarks form a color singlet background, and excitations are 
deviations from it. Quarks and quark- holes always appear together as confined 
excitations. 

To what extent can we generalize these pictures to the higher dimensions? 
The main changes can be found in the phase space allowed for colored qq 
fiuctuations which screen the exchange of gluons. 

In (1+1) dimensions, phase space near the edge of occupied states are always 
the same in the Dirac and Fermi sea. Thus the strength always looks same in 
a whole density region. This explains why an excitation energy near the Fermi 
surface does not change even at asymptotically high density. 

In spatial dimensions larger than one, the phase space for fiuctuation modes 
increases as density does. Eventually interactions between colors are strongly 
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screened, or squeezed color fluxes dissociate completely. Then weak coupling 
methods become enough to describe not only bulk quantities but also excita- 
tion modes. 

Currently we have only a rough estimate of such density for the (3+1) di- 
mensional system, jj, ~ A'^^.^/^Aqcd, paramerically larger than density where 
baryons overlap, /i ~ Aqcd- The relevance of the present discussions to (3+1) 
dimensions depends on how widely these two scales are separated after in- 
cluding all other numerical factors. Attempts to estimate such a width can be 
found in several model studies [30,31,32]. 

So far we have considered only high density regime in which the concept of 
Quarkyonic matter is most clearly illustrated. Now let us consider how such 
high density regime is connected to lower density regime where the presence 
of nuclear matter is important. 

Perhaps it is good to start with seeing how the picture of the baryonic Fermi 
sea in nuclear matter breaks down as density increases. The main problem is 
the difficulty to maintain the quasi-particle picture of baryons. Typical interac- 
tions among baryons get stronger at shorter distance, or in harder momentum 
transfer processes. This means that baryons near the Fermi surface and deeply 
inside of the Fermi sea can strongly affect one another. This badly destroys 
the concept of the baryonic Fermi sea at higher densitjQ 

After all, we have no right to stick to the Fock space made of conventional 
baryons made of quarks. The Fock space including baryonic objects made 
of 2Ai"c, 3A'"c, ■ ■ ■ quarks, which are also color singlet, may become equally im- 
portant. This is so, because energetically such states may be easily reached 
by strong forces among ordinary baryons as well as by large collision rates at 
high baryon density. 

This signals that effective degrees of freedom should be switched from baryons 
to quarks, in order to temper the growth of residual forces at higher density. At 
shorter distance, most of quarks feel weaker interactions, so they are reasonable 
alternatives forming the Fermi sea. Hence we expect that the transition from 
nuclear to Quarkyonic matter smoothly proceeds inside of the Fermi sea, while 
in both matters excitations remain confined. 

Here we should stress again that the basic ingredients of the above arguments 
are properties of interactions - whether composite objects overlap or not, is 
the secondary issue. For instance, if the baryon based descriptions gave smaller 



These arguments do not depend on whether is odd or even, or whether 
baryons are composite fermions or bosons. In high density nuclear matter, we start 
to observe their internal structures made of fermions. Needless to say, such quantum 
statics are important in dilute nuclear matter. 
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corrections of interactions than those in quark descriptions, we could continue 
to use the picture of the baryonic Fermi sea even at very high density. Of course 
such a situation is very unhkely in QCD. Needless to say, knowledges about 
baryon-baryon interactions at high density are crucial for further arguments. 

In summary, Quarkyonic matter differs from conventional quark matter in the 
excitation modes, and should be distinguished from nuclear matter by bulk 
quantities. 



3 The non-Abelian bosonization of QCD2 (One flavor) 



From this section, we start computations using the non-Abelian bosonization. 
The utility of the bosonized form is that we can compute several quantities in 
terms of quark number and color densities, emphasizing how a quark chemical 
potential acts very differently on quark numbers and colors. Another utility 
is that we can express explicit relations between chiral and quark number 
densities which are useful to consider issues of the chiral symmetry break- 
ing/restoration in quark matter. 



3. 1 Preliminaries 



We review the non-Abelian bosonization rules, to the extent which is necessary 
for our arguments (For quick introduction, [24]). The rules take compact forms 
by using lightcone coordinates x± = /2 = xq^Xi, and d± = 25^ = do±di. 
The currents for U{1) and SU{Nc) chiral charges are 



47r 47r 

= ^ii[hd.hHA] = -.i^ltA^-: , 4 = T^-irih^d+htA] = ■.i^itA^P+: , (8) 

ZTT ZTT 

with fields U {x) and h{x) which are frequently written as U {x) — e"^^^) and 
h{x) — e''^^'^^^^\ We have normal ordered currents and subtracted infinite 
constants. These currents are related to J^^ and J^5 as 






J.5 = e,,r =\ - / . (9) 



In chiral limit, they satisfy d^J^ = d'^J^^ = d±J^ = 0. Similar relations hold 
for J"^. The operator responsible for the chiral density is 
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r+^t = c(M) (/i«'')m(C/)m (10) 

where a, b is color indices of fermions. Here c(M) is a renormalization constant, 
and a subscript M means that operators are normal ordered at a scale M [33] . 
Note that the chiral density and the quark number density are characterized 
by common fields U and h, so once we know behaviors of these fields, we also 
know both densities simultaneously. 

The kinetic part of the fermion becomes 



/ d'x i;{x)zm^) S]:i'l[U] + SZfih], (11) 

where is an action of free massless boson, 

Sk^lP] = ^J d'x{d,U^d^U), (12) 
and S^J^[h] is Wess-Zumino-Novikov-Witten (WZNW) action [21], 



Sr^'il] - ^tr 



j (fx dlWl^'^ j (fx e'"'^{l^d^l){l^dj){l^dxl) 



(13) 



with level k — 1, and / = h. While the /3 function of the non-linear a model 
shows an asymptotic free behavior, the WZW term cancels out it. Thus the 
action is conformal at quantum level. A scale will be introduced by confining 
current-current interactions. 



Now we found one of the utilities of the non-Abclian bosonization. The colored 
currents include only colored boson fields h{x), and do not depend on other 
denisties in chiral limit: The actions for U (1) quark number and color charge 
sectors decouple. Only the color sector includes the dimensionful quantity, i.e., 
the gauge coupling constant with dimension one. The decoupled U{1) action 
is described by the free bosons and is conformal. 



3.2 The infrared fluctuations and the spontaneous symmetry breaking 



Before proceeding to our main discussions, first let us specify peculiarities in 
(1+1) dimensions. We will shortly take a glance at the relationship between the 
conformal property in the U{1) action and Coleman-Marmin- Wagner theorem: 
the spontaneous symmetry breaking is forbidden in (1+1) dimensions [14]. 
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Because of the separation of the U{1) and color sectors, the chiral condensate 
can be written in the factorized form, 



(VS+V-) = c{M){tT[h]M){UM). (14) 

This matrix element vanishes due to the infrared divergence in the propagator 
of the massless C/(l) boson which characterize the phase degrees of freedom in 
the chiral space. Intuitively, this reflects that phase fields in the chiral space 
rotate rapidly, without taking any particular direction. An explicit expression 

is (/3 = (47r/iVe)i/2) [33] 



Witten first observed this sort of A^c dependence in (1+1) dimensional corre- 
lation functions of chiral operators, which behave as ~ l/\x — y\^^^'' [13]. Here 
A is a propagator of the boson ip, and is the mass of (f which should be 
taken to be zero at the end of the calculations. 

This expression is one example which illustrates subtleties in the limiting 
order, large Ac and chiral limit. If we start with the strict chiral limit, we 
must first take limit before taking the large Ac limit. Thus the 

matrix element vanishes. 

The exception is the case with finite volume or with the infrared momentum 
cutoff. This is a typical situation studied in seminal works. They essentially 
yield similar effects to the introduction of m^. Once we regularize these in- 
frared fiuctuations, we can simply take the large A^c limit. 



2 2 / \ i Nc 



(^+^_) = c{M){tT[h]M) X {-^\ c{M){tr[h]M). (16) 

Here the dependence disappears as A^c — >^ oo- The renormalization constant 
and the value of (tr[/i]M) will be fixed in the next subsection. In the following 
we will frequently omit a subscript M if it does not play an essential role. 

Finally we would like to recall Witten's argument [13] that phase fields in 
(1-1-1) dimensions should not be identified as Goldstone bosons. While both of 
them share common features as phase fiuctuations, they play quite different 
roles in the chiral limit. The former belongs to the ground state properties, 
while the latter appears as an excitation from the vacuum. This aspect be- 
comes a little bit vague once we introduced an explicit breaking of the chiral 
symmetry, since phase fiuctuations become physical excitation modes after 
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the vacuum chose a particular direction in the chiral space. In what follows, 
we wiU not argue issues which arc very sensitive to this conceptual issues. We 
will regard phase fluctuations as excitation modes. 

3.3 A mass perturbation 

Now let us introduce a current quark mass term 

jC^ = -mg{ip+il)_ + = -TTiq X c{M){tT[h]U + tT[h^]U^). (17) 

which exphcitly breaks the conformal symmetry in the U{1) sector. Besides 
the infrared regularization, the mass term couples the U{1) and color sectors, 
thus it is no longer possible to investigate these sectors separately. Thus from 
now we rely on the large Nc limit. Then we can apply a probe approximation 
which significantly simplifies the arguments [8]. 

The point is that while the U{1) sector provides 0{Nc) contributions at most, 
those from quantum fluctuations in the color sector is 0{N^) because of 0{N^) 
degrees of freedom. The confining interactions generate a mass gap for the 
colored particles, so we can replace trh with (tih) in low energy. The value of 
(trh) is well-approximated by that computed without a current quark mass 
which is expected to act as a perturbation to the mass gap ~ Aqcd- 

We analyze the C/(l) sector under such a background. The effective action for 
low energy excitations is 



where we have subtracted a const ant to normalize the action. We have also 
used (trh) — (trh^) and defined = —rUg x 87rc(M)(tr[/i])/A^c for later 
convenience. Of course we can arrive at the similar Lagrangian following usual 
steps of the bottom up construction with the chiral symmetry. Only difference 
is that the above treatment clarifies roles of the colored sector. 

Now we take the expression U = e'-^^'^^^"'^^^'^'^ to canonically normalize the 
kinetic term for quantum excitations. Then the action for (p becomes 




(18) 




(19) 
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At large A'^c, should coincides with the mass of the lowest mode of the 
Bethe-Salpeter equations, computed by 't Hooft. He computed the pole mass 
of the channel J_ in the lightcone gauge. In a bosonized form, the long range 
behavior of the corresponding correlator is {x (g'^Nc)^^^'^) 



{j-{^)Am ^ {dM^)d-m), (20) 

so the pole should coincides with m^. Matching 't Hooft 's result with that in 
the bosonized form. 



= = (21) 



'C 



then we can fix c(M)(tr[/i]). Using these constants, we can see the value of the 
chiral condensate. Since regulates infrared behaviors of ip, the U{1) sector 
now has a nonvanishing expectation value. Thus we have 



{^l^^^;)^2c{M){tr[h])^-N^f-^. (22) 

This result was originally derived by Zhitnitsky using the current algebra plus 
the operator product expansion [9]. We reproduce this result simply because 
the bosonization procedures correctly implement the approximate chiral sym- 
metry. For later convenience, we write a pion mass in terms of the chiral 
condensate 



(23) 



4 One baryon, two baryons, and finite baryon density (One fiavor) 

4-1 One baryon as a topological object 



With a boundary condition for a topological charge, we can find a solution for 
coherent configurations U = e**^ regarding (f) as 0{N^) quantity. The action 
becomes the sine-Gordon model. 



^^^[0] = d^x[(9^0)^ + 2mj(cos0-l)], (24) 
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whose properties are well-investigated [7]. At large A^'c, stationary phase ap- 
proximation is applicable due to an overall A'";, factor of the action. Using 
Bogomol'nyi's trick, we have the lower bound of the energy, 



sm 



2t[ j0(-oo) 2 



>|A^b|x — ^ (25) 

where 0(±oo) must be 27rx integer. The equality holds only when the first 
bracket becomes zero. It is satisfied only by the topological charge one solution, 
whose configuration and energy are. 



(j){z; zo) = 4tan-i g-'^-^^-^"), E = (26) 

TT 

where zq is a free parameter to characterize a center of a single baryon. Here 
(j) goes to as z — )■ cxD and to 27r as 2; — ?■ — 00, as shown in the left panel of 
Fig.4. The baryon density is localized around z = zq a^s (the right panel of 
Fig.4) 



Jn 



Oz 

27r 



At: cosh (m^{z — Zq)) 
By integrating Jq with respect to z, we can see its quark number is A";,. 



(27) 



With expressions (26) and (27), the relationship between a quark number and 
chiral spirals is explicit. The scalar and pseudoscalar chiral densities are (here 
we attach a minus sign in front of A = |(V'V')vac| since ('(/''(/') vac < 0), 



(|)/2jt (m = m' ) 
(m = 2m' ) 




2 

1.5 

1 

0.5 





Jq (X 43i/N^) (m = m' ) 
(.51 = 2m' ) 




Fig. 4. The behaviors of 4> field devided by 27r (left panel) and quark number current 
jo devided by Nc/Air (right panel) as functions of z with some arbitrary dimensionful 
unit, m'~^. We chose zq to be 0. For each plot, we took niip = m! and m^p = 2m' to 
see the mass dependence. 
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A cos0(z; 2:0), 
= —A sin 2:0) 



(28) 
(29) 



so we find one chiral spiral around a single baryon. It is natural to expect 
more chiral spirals when we have a larger number of baryons. We will see this 
explicitly later. 

For future references, perhaps it is useful to rephrase the emergence of chiral 
spirals in a different way, without emphasizing (1+1) dimensional operator 
relations. Note that (ip+ip-) expresses an average density of pairs co- 

moving in the —z {+z) direction. Baryons make different potentials for pairs 
moving in opposite directions. Accordingly there is a mismatch in densities 
(ip+ip-) and near the baryon, so that (ipi'j^ip) does not vanish. We 

can repeat similar arguments in higher dimensions, once we replace 75 with 
7o7j, whose eigenvalues characterize moving directions of particles in a similar 
way as 75 does in (1+1) dimensions [3,16]. 

Finally we would like to examine qualitative differences in solitonic construc- 
tions of baryons in (1+1) and (3+1) dimensions. In the former, computations 
can be closed within dynamics around the scale, <^ Aqcd- This is in sharp 
contrast to (3+1) dimensional cases [34]: To get stable solitons, we need to 
equate the leading order chiral Lagrangian with higher order of derivative 
terms — they contribute equally, so that we do not know how to truncate 
expansions [35] . Another difference is that the former can be constructed from 
U{1) bosons only, while the latter is made of flavored pions. 

Presumably the simplicity in (1+1) dimensional construction comes from the 
fact that we see remnants of the charge-color separation. Such a separation 
is broken only via a term proportional to a small current quark mass, and 
the dynamics of colors only plays a very indirect role. Thus we expect that 
baryons discussed here have little to do with the picture of A^^c-quarks bound 
by the color fluxes. 



Fig. 5. The behaviors of scalar (S) {ipip) and pseudo scalar (PS) {ipi'j5ip) chiral 
condensates as functions of z. The unit of condensates is A. Left panel: 3D plot of 
the chiral spiral. Right panel: Separate plots of condensates. 
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4-2 Two baryons and interactions 



As we already saw in the one baryon sector, Bogomol'nyi's bound is satisfied 
only by single baryon solution. It means that the energy of two baryons is 
always larger than twice of single baryon energy, 



Eb=2 > 2Eb=i. (30) 

The additional energetic cost comes from the repulsive interaction between 
two baryons. Its asymptotic form is {R is the distance between two baryons) 



V{R) ~ N^m^e-'^^^. (31) 

When R ~ the strength of repulsion becomes the same order as the 

single baryon mass. 

How generic is this sort of the strong repulsive force in general dimensions? In 
QCD4, we know that the u meson exchanges in nuclear forces arc strong and 
repulsive. One interpretation of the strength is that it comes from the quark 
number of 0{Nc) which produces a large number of U{1) mesons of 0{Nc). 
The repulsive feature presumably comes from an additive property of U{1) 
charges, in sharp contrast to other non-Abelian charges such as isospins. 

We expect that this viewpoint be increasingly important at higher quark den- 
sity. While non-Abelian charge densities such as colors or isospins cancel one 
another and remain at 0(1), the enhancement of the quark number density 
is not stopped. Then strong repulsive interactions become increasingly impor- 
tant, and then completely change baryon structures. Eventually computations 
based on the quasi-particle picture of baryons break down, implying that the 
quarks become reasonable effective degrees of freedom. 



4.3 Finite density (one flavor) 



Since two baryon interactions are repulsive, baryonic matter starts to appear 
from the critical quark chemical potential, Hc = Mb/N^ = 2m^/'n ~ 0.7m^. 
In the following, the canonical approach is used. Since we have a large number 
of quarks, the leading contributions can be computed by studying coherent 
field configurations. With the coherent field expression U ~ e*"^, a total quark 
number of the system per period L is 



- I dzJQ[(f)] ^ - 
Jo 



, dz dzd — 

2ti Jo 2ti 



<p{L) - m 



(32) 
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Below we take (piL) — 0(0) = —2tiNb so that hb = Nb/L represents a baryon 
densit}p] We will use pf = ttub in the following. 

With this constraint of a finite quark number, we first determine a static 
background in the large Nc stationary phase approximation. The Hamiltonian 
for static coherent fields ij^ 



E\ 




-V. 



^{2ml+p}) + dzlid,^)' - 2mjcos(0 - 2pfz)V (33) 



where we have decomposed (f){z) field into 4>o{z) = —2pfZ and (f>{z). So 
satisfies a periodic boundary condition, 0(0) = 4>{L), forming the manifold 5*1. 
To find a solution for coherent configurations, we have to solve the following 
equation, 



d'^J - ml sin(0 - 2pfz) = 0. (34) 

The equation has been investigated extensively. It is known that solitonic 
configurations are general solutions which interpolate well-separated baryonic 
soliton configurations at low density and uniform quark number distributions 
at high density (For intensive discussions, see [36]). 

We will not repeat detailed analyses here but just quote some results about 
high density to assure our qualitative discussions in Sec. 2. 

In chiral limit, 0(z) = is the solution which means that the baryon density is 
uniform. The colored sector and the Lagrangian for quantum fluctuations are 
exactly the same as the vacuum. The chiral condensate behaves as {ipzfip±) = 
^g±i</> _ ^g±2ipjz^ form the chiral spirals. The coefficient A equals to the 



vacuum value The classical energy density is 



We require quantum fluctuations (p to be normalizable modes, so that the ex- 
pression (32) saturates the average baryon density. 

^ We assume that does not contain any constant, since it can be always absorbed 
by shifting 2pfZ — ^ 2pj(z — zq). 

It might seem strange that we get chiral spirals without any energetic minimiza- 
tion with respect to the chiral condensate. Actually, the modulus of chiral conden- 
sates is served from the color sector, independently from properties of the U{1) 
sector. If one hopes to see energetic differences with or without chiral condensates, 
he must derive an expression of the energy in the color sector. 
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The density contributions to the energy density is just those of free fermions. 

Next let us consider the case of massive fermions in asymptotically high den- 
sity. If we rewrite Eq.(34) by scaling a variable as z' — '2'PfZ, 



d'f4> + (^^J sin(^' - 0) = 0, (36) 
so ~ {rriip/pfY, and we can organize expansions. An asymptotic behavior, 



2 



0^^)=^ (^)'sm(2p/^), (37) 

satisfies Eq.(34) up to 0{m'^/p'j). This expression illustrates that roles of mass 
become less relevant as density increases. 

With the expression (37), we can compute a number of quantities of interest 
(Summary of distributions can be found in Fig.6.). The oscillation of the chiral 
condensate is slightly deformed as 



{i/ji/j) ~ A cos((/)o + 0^^^) = A cos (^pfZ - ^ sm{2pfz)^ , (38) 

and {'ipij5ip) also oscillates in the similar way. The baryon density also acquires 
modulations. 



{Jo) ^ -^dAct>o + = ^ X 2P/|l - ^ cos(2p;^) |. (39) 

In average, this solution does not contribute to the total baryon charge, since 
it just fluctuates around zero. At high density, prominent structures in quark 
number density is eventually buried in the average quark number background. 
The average energy density for ^ = ^^^^ is 



_ Kp) 




+ 0(m>^)|, (40) 

where the second term comes from the vacuum constant, while the third term 
expresses corrections from the density wave modulations. 
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Quark number density Jq 
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1.5 
1 
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Pf/fic = 0.5 



Pf/Hc = 1-0 



''jO "'lO '/iV '/ix" '.ri 



O'', o"', o"'' 
''-ii-"' ''',>>"" 



Chiral condensate {ipip) 






Fig. 6. The behaviors of quark number density (upper panels) and chiral 
scalar density (low panels) as functions of z with different quark densities 
Pf/t^c = 0.5, 1.0, 1.5. (The units of z, jo, (qq) are in~^, Nc^d'^'^-, ^i respectively.) 



An easy way to interpret expressions (38), (39), and (40), is to employ a 
particle-hole picture in the quark Fermi sea. When baryons largely overlap, 
quark number distribution is almost uniform without inducing color density. 
Deviations from such distributions can be regarded as corrections caused by 
particle-hole degrees of freedom. Indeed, spatial modulations have a period 
l/2pf, reflecting that condensations are driven by co-moving particle- hole 
pairs near the Fermi surface. Note also that the shape of peak in the baryon 
density (39) is very different from that of the single baryon solution, implying 
that structual changes of each baryon are essential. 



Based on the above arguments, it seems better to regard baryonic config- 
urations just as quark number localizations and defects, made of coherent 
configurations of color white particle-hole pairs. The dynamics of baryonic 
configurations are described as coherent excitations of such pairs (this gives 
more general descriptions than baryon and baryon- hole excitations). One of 
the relevant excitations are kinetic motions. They are characterized by zero 
modes associated with the violation of the translational invariance. 



With combining all these pictures together, our system is naturally viewed as a 
system with the Fermi sea of weakly coupled quarks plus mesonic excitations 
of particle-hole type. These arguments can be consistently combined with 
qualitative discussions in Sec. 2 which are based on the fermionic language. 
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5 Two flavors 



In this section, we extend analyses to two flavor cases with u and d quarks, 
keeping the vector symmetry SU{2)v by taking equal current quark masses 
= f^u = f^d- While most of treatments are the same as before, it is inter- 
esting to see the possibility of the chiral spirals including the flavor rotations 
such as {'4'ij5T3ip) . This issue may have phenomenological relevance to stars 
via the spontaneous generation of magnetic fields of the QCD scale. 

According to Vafa-Witten theorem [37] at zero density, the vector-like theories 
with current quark masses can not have the flavor symmetry breaking. How- 
ever, several conditions used in that proof can not be applied in the presence 
of the chemical potential, so this problem is interesting in its own right. In the 
following, we will not turn on electromagnetic interactions, but just see what 
happens in the QCD sector. 

We will see that, in the present model, the chiral spirals mainly occur between 
and ('0i75'0). The rotations including {ipi'y5Taip) can also occur, but it 
is order of rri^/ compared to U{1) channel, so will eventually disappear at 
high density. 



5. 1 Vacuum 



When we consider the multi-flavors, we have only to add slight modifications 
to previous treatments. The coefficients of U{1) and color currents change as. 



J_(x) = i-^Ud^U\ J^{x) = ij^iY[hd^hhAl (41) 
and flavor currents are defined as 



jL{x) — i—^iAgd- 



^-1 



' ^^^^^ ^ ^2n^^ 



(42) 



where the g \s a, matrix field and tj is a flavor matrix, with normalization 
tr [t/T/'] = 25///. The operators for the chiral density are 



r^i^'i = c(M) (/i'^'')m(C/)m(5''')m, (43) 

where i, j, . . . are flavor indices for the fundamental representation. The action 
for the fermion kinetic term is 
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(44) 



The treatments of the colored interactions are the same as before, and in chiral 
hmit, there is the charge-flavor-color separation. 

Below we will consider Nf — 2, massive quarks. With current quark masses, 
we have the following terms after applying a probe approximation and using 
tr[^]=tr[^t], 



^ ^^tr[g]{U + ^t) ^ ^^2"^^ cos n cos 0, (45) 
where we parameterized U — e^'^, and for the matrix field. 



g = e^^f^f = cos n + iilfTf sin U (H^ = U}, flf = Uf/Il). (46) 
This simple expression holds for Nf — 2, since {t/,t//} — 25ff/. 

From the form of Eq.(45), we can see that the trivial vacuum is satisfied when 
(0, n) = (0, 0) and (tt, tt) (modulo 27r), which maximizes the absolute value of 
the chiral condensate. On the other hand, the flavor or parity violating chiral 
condensates are 



^i-f^ip ~ tr[^](e^'^-e-^'^) ~ cos H sin 0, (47) 
V^T/V' ~ (tr[^T/]e^'^ + tr[^V/]e-*'^) ~ ft/ sin H sin 0, (48) 
^pi-f^Tf^ ~ (tr[5r/]e^^ - tr[5V/]e-^'^) ~ ft/sinHcost/), (49) 

becomes zero. These condensates compete with the usual chiral condensate. 
The usual chiral condensate wins because of the mass term. 

The action for quantum excitations can be obtained by expanding flelds around 
this vacuum. Taking U = e^^^'^/^^^f)'^'^ and g = e^^^'^/^'^)'^''^/^/, we have 



S^%,nf]^Jd'xl{{d,^r + {d,7rfr-ml{^' + 7r}))+0{l/^^^^^^ (50) 

The masses of bosons in U{1) and SU{2) are degenerate as a consequence 
of large N^.. At finite N^, there are higher order interaction terms in SU{2) 
sector, and they destroys such a degeneracy. 
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5.2 Finite density 



Calculations of finite density case can be done as before. We again consider the 
canonical approach, and take the average isospin density zero. Then we have 
a periodic boundary condition, 11/(0) — 11/ (L). The static energy functional 
for coherent field configurations is 

E[^,Uf] = AT.ATf ^ + Ekin + (51) 
where (below we will explicitly substitute A^f — 2) 

EU4>,'^f] = ^ J dz{{d,^f + {d,nf + {d^Uff sin' n}. (52) 

In chiral limit, it is clear that E^i^ can be minimized by simply taking = 
n = 0. Thus (j) — 2pfZ, so chiral spirals appear only in the C/(l) sector. This 
reflects that an external source (a quark number constraint) is put only in the 
U{1) sector. 

Now let us see what happens if we include mass terms without exphcitly 
breaking SU{2)v- The energy density from the mass term is 



Ejn[(l),Tif] — — ^ dz cos{2pfZ + (j)) cos U 

Stt j 



167r 



2 

^ / dz(^cos{2pfZ + + n) + cos(2p/^ + - H)) . (53) 



We will assume that and 11 do not contain any constant since they can 
be absorbed by shifting coordinate z. And it is more enlightening to rewrite 
kinetic terms. 



TV TV r 2 2 

Ekin[^,Uf] = -^J dz{{dS + U)) +(a,(0-n)) +(d,flf)'sm'u}(5A) 

From this form, dztlf = reduces the energy, so we choose IT = lis everywhere. 
Then we have two decoupled sine-Gordon models with finite density, whose 
variables are ± lis. We have only to borrow results of 1-flavor sine-Gordon 
model at finite density. 

As in 1-flavor case, amplitudes of and lis are ~ m^//i^ at high density. 
From expressions in Eqs.(45) and (47)-(49), we found the following asymptotic 
behavior of condensates. 
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(V'V'}, (V'«7^'0) « Acosn ~ A X (l + 0(m^///)2), (55) 
i'lpTfijj), {'ipi-f^Tfijj) oc Asinll ~ Ax 0{m^/iif, (56) 



so the chiral spirals mainly occur in the U{1) sector, while spirals with flavor 
breakings eventually disappear at high density. On the other hand, in relatively 
low density, we have solitonic lattices with mixture of a quark number and 
isospins. 



6 Terms disturbing chiral spirals 

In the preceding sections, we saw that the chiral symmetry is broken by chiral 
spirals when we have mechanisms to keep its modulus nonzero. We also ob- 
served that an explicit chiral symmetry breaking, i.e., a mass term, disturbs 
the chiral spirals, and such effects eventually disappear at high density. The 
purpose of this section is to generalize the above observation in such a way 
that we can classify several chiral effective models. 

First let us rephrase the role of a mass in QCD2 in the fermionic language, 

-mqjdz + (57) 

Since we have the Fermi sea, it is more natural to measure momenta and 
energies from the Fermi surface. Replacing 'il>± = tjj'^e^'^'^^ , we have 

-rriqjdz (VS^VVe^'"^ + i^'+i^'-e-^'"^) . (58) 

Suppose that the ground state gives (V'^pV'i) — const. This corresponds to chi- 
ral spiral solutions or uniform baryon distributions. For such a state, oscillating 
factors 6=*=^*^^ wash out energetic contributions from Eq.(58). This is what we 
actually observed in Eqs.(39) and (40). If ruq is comparable to /x, the above 
state is not the ground state and we should search solitonic configurations. 
Thus the mass term disturbs the formation of the chiral spirals. 

These arguments are useful to interpret results of other chiral models. For 
instance, 4-Fermi interaction terms in the chiral Gross-Neveu model (NJL2) 
with the continuous chiral symmetry are 

Hint j dz {{iP^f + (Vii75^)') ^^G j dz ((Vi+V'-)(VS-V'+))- (59) 

This form is unchanged when we rewrite fermion fields as 'ip± = ip'^e^"^^^ 
since oscillating factors cancel. Therefore the model at finite density takes 
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the same form as that in vacuum, except the constant energy shift associated 
with changes of the fcrmion fields. The chiral spirals or uniform quark number 
distributions start to appear for an infinitesimal chemical potential - there is 
no phase transition with increasing density. 

On the other hand, the discrete Gross-Neveu model takes the following inter- 
action term. 

The second and third terms have oscillating factors, and disturbs the formation 
of the chiral spirals. Thus chiral spirals and baryons can appear only after 
exceeds some critical value - there is a phase transition in contrast to the 
continuous model. At larger density, oscillating terms become less relevant, 
and the discrete model approaches to the continuous one (except the strength 
of coupling constant). Thus at high density, this model can be regarded as 
NJL2. 

These aspects are sometimes not manifest in the mean field energy functionals. 
Let us consider the discrete Gross-Neveu model. Typically one takes ansatz 
such as ~ 2A cos{qz) = S. This can be interpretted as 

{^^f = ((V;_V;+ - Ae*^") + (V^+V^_ - Ae-*^^) + S)\ (61) 

But the confusing point is that the final expression of the mean-field energy 
functional is characterized by {ipip) only, and roles of {i/ji'-f^ilj) condensate and 
chiral spiral structures are hidden. To see it, we have to explicitly calculate the 
condensate, or to construct the effective potential by inserting an infinitesimal 
external field. 

Finally let us apply these arguments to the chiral spirals or crystal struc- 
tures in higher dimensional systems. Consider states near the region (p^ ~ 
±/x, pt ~ Ot). By redefining fields * e^'*^'^"^'* = e^^^^^^i (^ is a fermion 
field in higher dimensions) , the fermionic kinetic term with a chemical poten- 
tial becomes 

£fcm ^ ^^id_\i,^ + + *+i^T*-e~^*''^ + *_^^T*+e^'''^ (62) 

so transverse kintetic terms disturbs the chiral spirals between (^'^') and 
(^7° 7^^). This field redefinition is useful as far as we consider only modes 
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with pt ^ 2/i. It means that condensations can occur only in the hmited 
region of the transverse momentum. 

Besides the aforementioned roles, transverse kinetic terms include coupling 
between and fields, and break the charge-color separation. This is 
consistent with the fact that in dimensions larger than one, increase of quark 
density affects the color sector through screening effects. 



7 Summary 

In this paper, we have utilized QCD2 to illustrate concepts of Quarkyonic 
matter. While Quarkyonic matter should differ from conventional quark mat- 
ter in the excitation modes, it should be distinguished from nuclear matter 
by bulk quantities. These aspects can be seen in quark matter in QCD2, as 
a consequence that quark chemical potential acts very differently on a quark 
number and color density. 

The nonperturbative dynamics near the Fermi surface play deterministic roles 
to classify the phase structure, since bulk contributions computed by weak 
coupling methods are approximately common for different phases. Whether 
excitations are deconfined quarks and gluons, or confined hadrons and glue- 
balls, is a key issue to understand dynamical phenomena in cold quark matter. 

Another relevant topic discussed in this paper was the inhomogeneous distri- 
butions of the chiral condensates and quark number densities. In our opinion, 
such inhomogeneous descriptions have potential relevance since they might 
smoothly interpolate the Fermi surface of Quarkyonic matter and very dense 
nuclear mattei p^ In both phases, excitations are confined, and chiral symme- 
try is broken. 

Several other important effects have not been addressed. In particular, we have 
not discussed the color superconductivity [38] for which a number of colors 

To avoid confusions, we have to emphasize that this crystalization is very dif- 
ferent from that occured in the Skyrme crystals [39,40] or its holographic QCD 
version [41]. The latter appears at very low density, because of long range attractive 
potential of 0{Nc) much larger than kinetic energy of P^/2M7v ~ l/iVc (Ref. [42] 
conjectured that this problem arises from the overestimated nucleon axial charge 
QA in the counting). This description contradicts with the real nuclear matter 
at low density, which are liquid-like rather than solid-like. In contrast, here we are 
discussing high density where repulsive hard cores of nucleons of 0{Nc) start to 
overlap, and nucleons can not move easily. This is the region where the nucleon 
Fermi sea starts to transform to the quark Fermi sea. The understanding of this 
region has not been reached. 
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are important. Since the formation of colored diquark condensates competes 
with the presence of confining forces, they must be taken into account simul- 
taneously for more realistic considerations than those given in this paper. 
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